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Abstract. The present knowledge on four-fermion production in electron-positron 
collisions is reviewed, with emphasis on W boson physics. Different methods to 
extract Mw from the data are presented and the role of QCD loop corrections 
discussed. 

1 Introduction 

Four-fermion processes represent the experimentally measured signal for W 
boson physics at LEP2 and NLC. In fact, the produced W's always decay, 
giving four-fermion final states [1]. One is primarily interested in measuring 
the W mass [2], but also measurements of the trilinear gauge boson cou- 
plings [3] and the four-fermion cross sections [4] provide useful information. 
An accurate determination of Mw can be combined with the precision data 
coming from LEP1 to further constrain the Standard Model of the Elec- 
troweak Interactions. In fact, a global fit to the LEP1 data [5] predicts the 
W mass with an error of the same order as the expected final LEP2 error (35 
MeV), while NLC will presumably reduce that error down to 15 MeV [6]. A 
strong discrepancy between predicted and measured Mw would be a signal 
for new physics. Alternatively, an improvement on the measurement of Mw 
can significantly tighten the present bounds on the Higgs mass through the 
relation 

Q = ^ 1 fl } 

" V2M W (1 - M W /M 2 Z ) 1 - A r (m t , m H ) ' K ' 

where A r is a calculable contribution coming from radiative corrections. 

Several tree-level four-fermion codes are available [4] . Electroweak radia- 
tive corrections are usually included at the leading log (LL) level, but very 
recently, non-factorizable QED corrections have been computed as well [7]. 
While tree level programs + LL corrections seem in general to be adequate 
to deal with LEP2 physics, further refinements, especially in the sector of the 
loop radiative corrections, are needed in view of the NLC precision physics. 

In this contribution, I focus my attention on two particular aspects of 
four-fermion W physics, namely Mw measurement and QCD corrections. 

In the next section, I describe two techniques to measure the W mass, 
and present a new method to extract Mw from the best measured variables. 
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When quarks are present in the final state, QCD loop contributions have to 
be included as well. Those corrections are discussed in the last section of the 
paper. 



2 M w measurement 

Two methods are mainly used to extract the W mass: the threshold method 
and the direct reconstruction technique [2]. In the first case the total W + W~ 
cross section is measured near threshold (161 GeV), where the sensitivity to 
Mw is stronger, and plotted as a function of the W mass. At LEP2, that 
gives M w = 80.40 ± 0.22 GeV [5]. 

The direct reconstruction method is applied at higher energy, where the 
statistics increases. It requires three steps: 

1. From the experimental data the invariant mass distribution is re- 
constructed. To improve the mass resolution, a constrained fit is usually 
performed event by event, assuming no Initial State Radiation (ISR) and 
equality between the invariant masses coming from different W's. 

2. A theoretical distribution is taken for -fe (usually a convolution of a 
Breit-Wigner with a Gaussian) and a mass M' w fitted. 

3. The value M' w is then corrected by Monte Carlo, from the bias introduced 
by the constrained fit, to get the reconstructed W mass M R with an error 
AM R . 

Such a measurement gave M w = 80.37 ± 0.19 GeV [5], in the LEP2 run at 
172 GeV. 

Recently, a new method has been proposed [8] (direct fit method), in 
which only the best measured quantities are used to extract the W mass. 
The idea is simple. Given a set of well measured quantities {<P} one computes, 
event by event, the theoretical probability Pi of getting the observed set of 
values {<I>i} for This is a function of M w and is given by the ratio of the 
differential cross section in those variables, divided by the total cross section 
in the experimental fiducial volume 

da 

P i {M w ) = ^-. (2) 



Given N observed events, the logarithm of the likelihood function L 

N N 

lna 

d<P. 



log L(M W ) = log HPi(M w ) = E lo gS-( M ^) -Nloga(M w ) (3) 



is distributed, for large N, as a quadratic function of Mw- The previous 
equation is then computed for different values of Mw and a parabola fitted, 
from which the reconstructed W mass Mr is obtained with an error AMr. 
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In order to construct a tool for the evaluation of Pi(M\y) one has to choose 
the set {</>} of accurately measured variables. Although one can always con- 
sider more sets {</>}, the following choices seem reasonable in practice [9] for 
different four-fermion final states: 

1. Semilcptonic case: q\q^lv 

la {cj>} = {E e ,n e ,n qi ,n q2 } 

lb {<p} — {Eg, Qi, fi qi , fi q2 , Eh}, where Eh is the total energy of the jets. 

2. Purely hadronic case: 

{<f>} = {n qi fi q2 n q3 ft qi }. 

3. Purely leptonic case: l\V\iiVi 

{<j>} = {E tl ,n il ,E t2 ,n l2 } 

Since eight variables determine an event when no ISR is present, sets la and 
3 would require one or two integrations, cases lb and 2 none. Including ISR 
adds two integrations. When jets cannot be assigned to specific quarks a 
folding over the various possibilities should be included. 

As an example of the direct fit method, I show, in table 1, the recon- 
structed masses obtained by fitting a sample of 1600 EXCALIBUR [1] CC3 
events [4] including ISR, generated with an input mass of Mw = 80.23 GeV, 
at = 190 GeV. In all four cases, the W mass is correctly reconstructed. 



la 


lb 


2 




80.238 ± 0.049 


80.238 ± 0.032 


80.255 ± 0.036 


80.209 ± 0.077 



Table 1. Reconstructed W mass (GeV) with four choices of the set {<&} (see text). 



Cases lb and 2 give better errors, because less information is integrated 
out. Conversely, in the leptonic case 3, the error is worse, since a large part 
of the kincmatical information is actually missing. 

In the direct reconstruction method, one wants to keep as much infor- 
mation as possible, also preferably photon momenta, in order to reconstruct 
the kinematics event by event. On the contrary, in the direct fit method, one 
has to integrate over all information that is not well determined, in partic- 
ular ISR. Because of the fact that the integral over the px distribution of 
ISR photons is theoretically better known than the distribution itself, one 
expects the details of the radiation to matter less in the direct fit method. 
That can be viewed as an advantage with respect to the direct reconstruction 
technique. 

However, a last remark is in order. All numbers presented in table 1 refer 
to the partonic level, without inclusion of hadronization effects and detector 
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resolution. Therefore one still has to prove that the fitting procedure survives 
those effects. This question is currently under investigation [10]. 

It is also clear that the direct fit method is not only applicable to measure 
Mw, but can be used, in principle, to extract any parameter - as IV or a 
set of anomalous trilinear gauge couplings (TGCs) - from the data sample. 
Not surprisingly, the whole strategy for the direct fit method has been first 
discussed in the final report of the Workshop on Physics at LEP2, in the 
context of TGCs determination [3] . 

3 QCD corrections 

QCD loop corrections to four-fermion production in e + e~ collisions can be 
divided in two classes, namely QCD corrections to 0(a 2 a 2 ) and 0(a 4 ) pro- 
cesses, respectively. 

The first corrections appear as 0{a 2 s ) contributions to four-jet produc- 
tion via QCD [11], while the second ones are relevant for studying W boson 
physics, and, more in general, semileptonic four-fermion processes and fully 
hadronic final states mediated by electroweak bosons. 

I shall concentrate here on the latter contributions, which can all be ob- 
tained by defining suitable combinations of loop diagrams plus real gluon 
radiation, as shown in ref. [12]. The calculation is simplified a lot by using 
the reduction procedure presented in ref. [13]. 

In table 2, cross sections computed with the program in ref. [12] are 
presented, for the semileptonic process e + e~ — > ^v^ud. 





Born 


NLO 


nQCD 


161 GeV 


.24962 ± .00002 


.24760 ± .00002 


.24790 ± .00002 


175 GeV 


.96006 ± .00007 


.94519 ± .00007 


.94613 ± .00007 


190 GeV 


1.184003 ± .00009 


1.16681 ± .00009 


1.16766 ± .00008 


500 GeV 


.46970 ± .00006 


.47109 ± .00007 


.46131 ± .00006 



Table 2. Cross sections in pb for e + e — > /i v^ud with canonical cuts [4] . 



The exact calculation (NLO) is compared with a "naive" approach to strong 
radiative corrections (nQCD), where the QCD contributions are simply in- 
cluded through the substitutions 
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r„-iv(i + §*), „-<,(i + £i). ,4, 

For the direct reconstruction of the IT mass, the quantity 

(AM) — J (y / s^+ — 2 Mvk) da is relevant [4], where M\y is the input 

mass in the program. One gets, with canonical cuts at y/s = 175 GeV: 

(AM} N lo = - 0.5585 ± 0.0002 GeV 
(Z\M) nQC D = - 0.5583 ± 0.0002 GeV . (5) 

Also, one can show that the angular distributions are distorted, in the NLO 
calculation, with respect to the nQCD prediction [12]. 

From the previous results it is clear that nQCD is adequate at LEP2 
for semileptonic processes, but exact calculations are important at NLC and 
for anomalous couplings studies, where the angular distributions matter to 
constrain the anomalous contributions. 

In table 3, I show results for the fully hadronic process e + e~ — > udsc at 
y/s = 175 GeV. The numbers are obtained by using the program described in 
ref. [14]. In case (a) only canonical cuts are applied. In (b) two reconstructed 
masses Mm and Mr2 are determined by minimizing the quantity A' M = 
(Mri - M w ) 2 + (Mr 2 - M w ) 2 , and a cut \M m - M w \ < 10 GeV is imposed. 
In (c) a smearing with a Gaussian with a 2 GeV width is introduced in 
addition, to mimic the experimental resolution. 



a(pb) 


(a) 


(b) 


(c) 


NLO 


1.1493(4) 


0.7895(5) 


0.7758(9) 


nQCD 


1.1069(3) 


1.0545(3) 


1.0479(3) 



Table 3. Cross sections for e + e — > udsc at \fs = 175 GeV. 



For the process at hand one gets, for case (b), 

(AM) NLO = - 0.2290 ± 0.0010 GeV 

(Z\M) nQC D = - 0.0635 ± 0.0004 GeV , (6) 

where now (AM) = ± J (M m + M m - 2 M w ) da. 

We are easily convinced, by the above results, that the naive QCD im- 
plementation fails in describing hadronic four-fermion final states at LEP2. 
In particular, the reduction in cross section (compare cases (a) and (b) in 
table 3) shows that many soft gluons are exchanged between decay products 
of different VF's that are not taken into account using nQCD. A failure of 
nQCD can also be proved at NLC energies [14]. 
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4 Conclusions 

Precise measurements of Mw are performed at LEP2, using the threshold 
method and the direct reconstruction of . 

An alternative technique is also available in which only well measured 
quantities are directly used to extract M w (and TGCs) from the data. 

QCD loop corrections to semileptonic 0(a A ) processes are well approxi- 
mated at LEP2, by nQCD, except for angular distributions, while the naive 
approach fails in describing hadronic four-fermion final states. 

All existing calculations have to be refined and radiative corrections better 
understood in view of the NLC precision measurements. 
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